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Abstract 

We describe conditions under which an appropriately-defined anisotropic 
Voronoi diagram of a set of sites in Euclidean space is guaranteed to be 
composed of connected cells in any number of dimensions. These condi- 
tions are natural for problems in optimization and approximation, and 
algorithms already exist to produce sets of sites that satisfy them. 

1 Introduction 

The anisotropic Voronoi diagram (AVD) is a fundamental data structure with 
wide practical application. In the definition of [9 , an AVD over a Riemannian 
manifold is the Voronoi diagram of a set of sites with respect to the geodesic 
distance associated with a Riemannian metric. From a formal viewpoint, this 
definition has several strengths. For instance, a simple argument can be used 
to show that the Voronoi regions of such an AVD are always connected. 

For problems in optimization, such as vector quantization on a Riemannian 
manifold, the above property makes it simpler to compute the AVD: since each 
Voronoi region is connected we do not have to search for possible disconnected 
(orphan) pieces elsewhere. For problems in approximation, where we are often 
more interested in the dual simplicial complex, the absence of orphan regions in 
the AVD makes it possible for strict approximation error bounds to be enforced 
on the dual and, in some cases, it can make it possible to characterize the 
asymptotic size of the approximation [71 [3] . Geodesic distance on a Riemannian 
manifold is, however, very expensive to compute, since it involves finding the 
shortest geodesic path between two points, and a practical AVD algorithm based 
on this distance is not currently feasible. 

By choosing a particular parametrization of a Riemannian manifold over a 
subset of Euclidean space, two approximations to the geodesic distance between 
two points naturally arise, both of which have been considered as basis for 
constructing AVDs. These approximations simplify the problem by considering 
the metric constant along any path between the two points, but use a different 
choice of constant metric along the path. In particular, to measure the distance 
between a given site and any point in the domain, one approximation evaluates 
the metric only at the site, while the other uses the metric at the point. They are 



described in [5] and [5]. We refer to them as Labelle/Shewchuk, and Du/Wang 
diagrams, respectively, or LS and DW diagrams for short. Although they are 
conceptually similar, guarantees of well-behaved-ness have only been shown only 
for the case of LS diagrams, and only in two dimensions [S]. These guarantees 
come in the form of a condition which, if satisfied, guarantees connectedness 
of Voronoi cells as well as the well-behavedness of its dual 2D triangulation 
(absence of triangle inversions). They also provide an iterative site-insertion 
algorithm that is guaranteed to produce a well-behaved output if run for long 
enough. 

In this paper we present conditions under which these anisotropic Voronoi 
diagrams are well-behaved. These conditions are simple and intuitive, hold for 
any number of dimensions, and apply to both LS and DW diagrams. 

The condition we describe requires that the generating set of sites form a 
sufficiently dense e-net (with e being sufficiently small) . This condition is quite 
natural, since it requires that the distribution of generating sites be roughly 
"uniform" with respect to the input metric. An iterative, greedy site-insertion 
algorithm [5] exists to compute sets that satisfy it, which works by iteratively 
inserting a new site at the point in the domain that is furthest from the cur- 
rent set of sites. Note that the net condition is natural both for optimization 
problems, such as quantization or clustering, where we are interested in the 
primal Voronoi diagram [3 Sill], as well as for approximation problems, where 
we are mostly interested in the dual simplicial complex. In the former case, 
optimal quantization sites have been shown to satisfy the net property [71 13] , 
while in the latter, nets have been used to construct asymptotically-optimal 
approximations [3 . 

The net condition is both a condition on the density of sites (cover property), 
and its relative distribution (packing property). While the packing property is 
clearly insufficient to guarantee that an AVD is well-behaved, we show that 
the density property alone can be insufficient in cases when the net condition 
is. That is, even though there may be sufficiently dense covers that produce 
well-behaved AVDs, there are e-nets that produce well-behaved AVDs where 
the e-cover property alone would not (the converse is not possible since every 
net is a cover). For sufficiently fine densities of sites, the combination of cover 
(density), and packing (relative position) properties is enough to ensure that 
both the DW and LS diagrams are well-behaved, in any number of dimensions. 
Bounds on the minimum density required are given which, perhaps surprisingly, 
do not depend on the dimension. 

2 Preliminary definitions 

Given an Euclidean domain C M", an ordinary Voronoi diagram of a set 
y C is a partition of the domain into regions whose points are closest to the 
same element in V. In the case of LS or DW diagrams, the function used to 
measure closeness is not the natural distance. If we assume that we are given 
a Riemannian metric over 0, with coordinates Q : — M"^" (where at each 
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p € 0, Qp is symmetric, positive definite), we define the functions: Dq {a,b) = 
[{a - by Qa{a - b)]^^'^ , and Dq"^ {a, b) = [{a - bYQb{a - b)]^^"^ {noteDQ{a,b) = 

DW 

Dq (6, a)). An LS diagram (resp. DW diagram) of a set V is the Voronoi di- 
agram of V with respect to the function D^ (resp. Dq" ). The associated 
Voronoi regions of a site v G V are, respectively, 

={pe 0\yw e V,Dq{v,p) < DQiw,p)} 

and 

= e 0|Vw e y,i?Q"(i;,p) < Dq\w,p)} 

Note that, because neither Dq nor Z3q are symmetric, we must follow some 
convention on the order of its arguments. In particular, we place an element of 
V always as first argument of Dq and Dq . We follow this argument order 
convention for the rest of the paper. 

Asymmetric £-net. In the remainder of the paper, it will be useful to consider 
e-nets with respect to a function that is not symmetric (e.g. Dq and Dq ). 
Since the original definition of an e-net assumes the use of a symmetric distance, 
we must slightly modify it in this asymmetric case. As before, we follow the 

DW 

convention of placing elements of the net always as first arguments to Dq or 




Definition 1. An asymmetric e-net with respect to a function D : x — > K 
is a set V Cz that satisfies: 

1. Vp G 0, mmD{v,p) < e. (asymmetric e-cover property) 

2. yv,w e V, D{v,w) > e or D{w,v) > e. (asymmetric e-packing property) 

These properties are analogous to those of a regular net, but not identi- 
cal (e.g. the packing property is weaker). Note that, even if D above is not 
symmetric, we can still compute an asymmetric net using the greedy algorithm 
of [6j, by being careful to follow the stated argument-order convention. A simple 
induction argument reveals that, in the asymmetric case, the algorithm of [B] 
always terminates by outputting a discrete set V that satisfies the properties 
of an asymmetric net. (In particular, the output of W\ will not, in general, 
satisfy the stronger version of the asymmetric packing property: Vv, w & V , 
D{v,w) > e and D{w,v) > e). 



3 Setup 

Assume that we are given as input a metric over n-dimensional Euclidean space 
(in coordinates: a field of symmetric, positive definite (PD) matrices Q : — >■ 
M"^"). At every point in 0, there is an eigen-decomposition Q = RKR*^ . Let 
the symmetric PD square root matrix of Q he M = Rh}/^R^. As in |8|, we 
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use this square root matrix to analyze the Voronoi diagrams under consider- 
ation. In contrast to [S], we do not consider any square-root matrix M' such 
that M'*M' = Q (ah square roots can be written as UM, where U is unitary: 
/7*C/ = /), but rather concentrate on the unique square root that is also sym- 
metric and positive definite. This distinction is important. In particular, as the 
following lemma shows, the symmetric PD square root has the same differen- 
tiability class as Q, and, in particular, it is continuous wherever Q is. 

Lemma 1. Given Q : G — )■ R"^", symmetric, positive definite over a 
compact domain 0, then the unique symmetric, positive definite matrix M that 
satisfies M^M ^ Q is M & C'' if and only if Q ^ 

Proof. If M G then clearly Q = M*Af ^C^. To show that the converse is also 
true assume Q ^ and consider, at any point p e 0, the eigen-decomposition 
Q — RAR*, where A is a diagonal matrix with positive entries. We claim that 
M = i?A^/^i?*. Clearly, M is symmetric, positive definite, and M'^AI = Q. Any 
other positive definite matrix M' with M'*M' = Q \s related to M by a unitary 
matrix as M' — UM. But, since M is symmetric, M' can only be symmetric if 
U = I. Therefore AI = i?A^/^i?* is the only square root of Q that is symmetric 
PD. 

Just like the Taylor expansion of the (scalar) square root function ^/x = 
X]^o (i-2n){nrf(4") ~ 1)"' which couvcrgcs for ja: — 1| < 1, we can define the 

series S — J2^a (i-2i^Ki!)^(4') (Q ~ -^)'' "^hich, by d'Alembert's ratio test [TU], is 
absolutely convergent for symmetric PD matrices with spectral radius p{Q) < 2 
since the absolute value of the coefficients in the series is decreasing, and the 
norm of (Q — /)* is exponentially decreasing with i. We show that if p{Q) < 2, 
then the series converges to M, the unique symmetric PD square root. Clearly, 
iS* is symmetric and, given an eigenvector of Q with eigenvalue h, it is elSe^ = 

,poo (-l)"(2ri)! n i yi _ ;l/2 Vp^irp ^ — M 

The derivatives Sn^ = J2i=o '^^''^ (Q ~ -^)' of the partial sums are 

defined for j — 0, . . . , fc, and have the same convergence region as the original 
series S. Since Q is symmetric PD, p{Q) < 2 and is compact, Q — I has 
maximum spectral norm pmax = maxpgn p(Q — /)<!. Thus the slowest rate of 
convergence occurs at points p € fl where p{p) = Pmax, and therefore the Sn^, 
j = 0, . . . , k converge uniformly inside 0. Uniform convergence implies that the 
limit of derivatives is the derivative of the limit 10] : lim„_>.oo Sn'^ = d'^^S' = 
d'^^^M. Therefore M is at least in the same differentiability class as Q. 

The above applies to matrices Q with p{Q) < 2. In the general case, scale 
Q by l/maxpg0 p{Qp), apply the lemma, and rescale back the resulting M by 
maxpg0 □ 

Because Q is spatially- varying, given two points a,b € 0, it will in general be 
Dq (a, 6) =/= Dq (6, a), and likewise for Dq . The amount of asymmetry will 
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depend on how different Qa is from Qb. In particular, by the sub- multiplicative 
property of the spectral matrix norm, 

DW , , , DW , 

Dq {a,b)/DQ {b,a) 



and a similar argument shows that i'g" (a, h)/D'^^ {b, a) > pm{Mi,M~^), where 
p{A) is the spectral norm and Pm{A) is the smallest of the absolute values of the 
eigenvalues of A (that is: p{A) — sup^-^Q || Ar||/||r|| and Pm{A) — inf^^^io ll^^ll/lkll) 
Note that, for square matrices, such as the ones considered here, the spectral 
norm is the same as the induced operator norm. 

Given a point a € 0, it is possible to bound the amount of asymmetry in 
Dq (resp. Dq ) inside an appropriately-defined neighborhood of a. To this 
end, we introduce the following definition, which is applicable to continuous 
metrics, whose square root M is, by Lemma [l| also continuous. 

Definition 2. The maximum variation of a metric Q is the smallest constant 
tr such that for all a, 6 G 0, it is 

p {MbM-' ~ I) < a ■ \\Maia ~ b)\\ 

where M is the symmetric, positive definite square root of Q. 

Loosely speaking, cr is a Lipschitz-type bound on the rate of variation of M 
relative to itself. In the sequel, it will be assumed that a is finite. In particular, 
this will always be the case if is compact. 

We can use the above definitions to find bounds on the asymmetry in the 
associated function Dq (resp. Dq ). The following lemma shows that, if a 
point b € is inside a certain neighborhood of point a e 0, and a is suffi- 
ciently small, then p{Mi,M^^) can be bounded from above, and Pm{Mi,M~-^) 
from below, which implies that Dq (a, b) and Dq (6, a) (resp. Dq (a, b) and 
Dq (fo, a)) must be similar. 

Lemma 2. Given e > and Q E with maximum variation a , then for all 
a,b € with \\Ma{a — b)\\ < e it is 

l-ea< pmiMbM-^) < ||Mb(a - b)\\/\\Ma{a - b)\\ < p{MbM-^) < 1 + ea 

Proof. The upper bound follows from the sub-multiplicative property of the 
spectral norm and the definition of a: 

\\M,ia - b)\\/\\Ma{a - b)\\ - \\M,M-Hlaia - 6)||/||A4(a - < p{M,M-') 

- p{MbM-^ -/ + /)< 1 + p{MbM-^ - I) 
< I + \\Maia - b)\\a < 1 + ea 



Dr 



\b,a)/DQ {aM = UUa-b)\\/\\AUa-b)\\ 



\\M,M-Hla{a^b)\\ 1 
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If li are the eigenvalues of a symmetric matrix A, then 



Pm{A + T) 



= min 1^,; + 1| > min ||1| — 



i i 



> min 1 — I Z,j I = 1 — max =■ 1 — p{A) 



and therefore 



\\Mb[a-b)\\l\\Ma{a-b)\ 



> 



> 



\\MbM-^Maia-b)\\/\\Ma{a-b)\\ 
p^MbM-^ -/ + /)>!- p{MbM-^ 
1- ||A/a(a- 6)|1ct > 1-ea 



I) 



□ 



In the remainder of the paper, we estabhsh conditions for anisotropic Voronoi 
diagrams to be orphan-free. In particular, the set of sites considered will be 
asymmetric e-nets, where e must be sufficiently small in relation to the constant 
(T above, which depends on the input metric and, vaguely speaking, provides an 
upper bound on the rate of change of the metric. 

Note, however, that, in practice, it is not necessary to compute a to find 
a sufRciently small e that guarantees that Voronoi regions are well-behaved. 
Instead, it is possible to simply run the greedy algorithm of [5], which in our 
case outputs asymmetric e-nets whose e decreases with each iteration, until the 
resulting Voronoi diagram is orphan free. This is because, at each iteration, 
the algorithm of |^ must compute the closest site to each point in the domain 
- a task that is equivalent to computing the Voronoi diagram for the current 
set of sites. Therefore checking at each stage whether the current diagram is 
orphan-free can simply be a by-product of the asymmetric e-net computation 
algorithm. The proofs in this paper simply guarantee that there is a small- 
enough e for which the resulting asymmetric e-net produces an orphan-free 
Voronoi diagram, and thus that the above algorithm stops (a proof involves 
the fact that the Voronoi radius of a site can be made as small as desired by 
simply introducing more sites in the e-net, which follows from the definition of 
the DW/LS distances). The precise bounds in Theorems [l] and [5] may also serve 
to give some indication of how small e will need to be. 

4 Orphan-free anisotropic Voronoi diagrams 

This section shows that, given a continuous metric, the associated DW and 
LS diagrams of an asymmetric e-net, for sufficiently small e, are orphan-free. 
Theorems [T] and [2] state conditions under which this holds. In particular the 
criteria for e to be sufficiently small will be a certain relation between e and ct, 
the maximum variation of the metric, which does not depend on the dimension. 
(Specifically, ea < 0.09868 for DW diagrams, and ea < 0.0584 for LS diagrams.) 

Since the (asymmetric) e-net property is a combination of an e-packing and 
e-cover properties, it is natural to consider whether any of these properties is. 
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Figure 1: If two sites (v and w) are very close-together, a small perturbation 
in Dq can cause an orphan (top left), even if the sites are very dense with 

DW 

respect to Dq (top right). If sites are sufficiently spaced apart relative to each 
other (bottom left), orphan regions need sufficiently-large (relative) fluctuations 

DW 

in Dq to appear. In this case, placing the sites more densely-packed with 

DW 

respect to Dq reduces the relative fluctuation, which eliminates the orphans 
(bottom right). 



by itself, sufficient to guarantee that DW and LS diagrams are well-behaved. 
Clearly, the e-packing property cannot be sufficient since any two sites form 
an e-packing for any s smaller than the distance between them. The case of 
e-covers is more subtle: for a particular choice of metric Q, it is possible that a 
sufficiently small e' exists such that every e'-cover produces well-behaved DW 
and LS diagrams. When we consider all possible choices of Q, however, the 
required e' may be arbitrarily small relative to a and, unlike for e-nets, there 
is no constant c > such that, for all choices of Q, an e-cover with ea < c is 
always guaranteed to produce orphan-free DW and LS diagrams. 

To see this, consider the diagram of Fig. [T] In the right column, the set of 
sites (black dots) is dense enough over so that Q is roughly constant inside 
each Voronoi region. There is, however, some small variation in Q. We can 
place two sites v,w € V very close together (but not coinciding) such that, even 
for small e, a very small change in Q away from them causes an orphan region 
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to appear (top right). A point in the orphan region, near the interface between 
v^s and ui's Voronoi regions, "sees" both v and w as being at approximately 
the same distance, but a very small variation in Q has made the points in the 
orphan region be slightly closer to v, even though they are surrounded by points 
that are slightly closer to w. Although for a particular choice of metric Q there 
may be a sufficiently small choice of e for which all e-covers are guaranteed to 
produce well-behaved DW and LS diagrams, the fact that the variation in Q 
described above can be arbitrarily small means that the requirement on e may 
be arbitrarily strict, depending on the choice of Q. 

4.1 Orphan-free Du/Wang diagrams 

DW 

Given a continuous metric Q and its associated distance Dq , we show here 

DW 

that an asymmetric e-net with respect to Dq , with sufficiently small e, has 
Voronoi regions that are always star-shaped (with respect to the Euclidean 
distance) from their generating sites and, in particular, they are connected. We 
prove this by showing that, if a point g g belongs to an orphan region of 
some site w G V^, as in the diagram (Fig. [2]a), then the segment qv connecting 
them must also belong to the Voronoi region of v, contradicting the fact that 
the Voronoi region that contains q is disconnected from v. 

More specifically, since q is in an orphan region, the segment (jv must contain 
a point p' that is closer to a different site w G V {p' is in the interior of the 
Voronoi region of w 7^ w). However, we show that p' cannot be closer to any 
w V, reaching a contradiction. In conclusion, qv must belong to the Voronoi 
region of v, and so q cannot be in an orphan region. Additionally, this shows 
that every Voronoi region must be star-shaped with respect to its generating 
site. 

The details of the proof follow. In particular, we will see that, if q is in an 
orphan region of v Cz V, and Q is continuous, then by the intermediate value 
theorem, there must be two points c,c' € qv that are at equal distance to v and 
to some other w &V. We use the e-cover and £-packing properties of V to show 
that, for sufficiently small e, the existence of such c, c' leads to a contradiction. 

Assume q & Ry and p' e R^, , where R^ , R^ are the Voronoi regions 
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of V and w respectively. Then there must be a point c e qp" between them that 
belongs to both Voronoi regions, such that Dq {v, c) = Dq {w, c). Since V is 

an e-cover, this point must also satisfy Dq (w,c) = Dq {w,c) < e. 

Consider now the parametrized segment p{l) = c(l — I) + vl with / e [0, 1]. 
Letting m ^ (v + w)/2, define the function 

,1 r. 



DW 



=£"'2 r« Kp(0)-^q {v,p{l))\^e-'{w-vyQ^(i){m^p{l)) 

which is plotted in Fig. [2jb. Note that (j){p{l)) is continuous with respect to I 
by virtue of the fact that, by Lemma [ij it is M S C*^. 

Since c is equidistant to v, w, it is 0(c) = 0. For / > 0, becomes negative 
at p' , since p' is closer to w than to v, and then becomes positive at v (since 
V is closer to v than to w). Because it is 4>{p') < and 4>{v) > 0, and since it 
is (j) G C'^, then there must be an intermediate point c' € p'v with 4'{c') = 0- 
Finally, because (/)(c) = (j){c') = 0, it is 

= f!^(0(c)-0(c'))=e-2fc4(«;-«)*[Q,,(m-c')-Qc(m-c)] 



= £-2^^ [(c - d)'Q,(w -v)-r(w- vYiQc' - Qc)(m - c')] 
1 1 c — c' 1 1 

= e-2(c - vYQ,{w -v)+ e-'^p^iw - «)*(Q,, - Q,)(m - c') 

||c-c'|| 

If we define a = |e"^fe|(w - vYiQc' - Qc){m - c')|, and /3 = \e-^{c ■ 



vYQciw — then 



||c - c',, 

We can reach a contradiction by showing that (/)(c) ~ (/'(c') does not vanish. To 
do this, it suffices to bound a from above, and j3 from below in such a way 
that their difference is always positive. In particular, we will see that a can be 
made arbitrarily small by requiring V to be an e-cover of sufficiently small e. 
To bound j3 from below, on the other hand, it is not sufficient for V to form a 
sufficiently dense cover. /? is sensitive to both the density of sites in V , as well 
as their relative distribution. It is, however, possible to find a sufficiently-high 
lower bound of j3 by requiring V to be an asymmetric e-net. The asymmetric 
e-net condition is therefore sufficient to bound both /3 from below, and a from 
above, in such a way as to ensure that (/)(c) — 0(c') doesn't vanish, creating a 
contradiction and concluding the proof. The following two lemmas provide the 
relevant bounds for a and /3. Auxiliary lemmas from the Appendix are used in 
the proofs. 

Lemma 3. Given an asymmetric e-cover V , ifv,w G V are Voronoi-neighbors, 
and it is q € R.^ and c, c' and m as described above, then 

a = le-' ll"'";! [w - vYiQc' - Qc){m - c')| < 2 {eaf + Aea 
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Proof. Since c is in the Voronoi regions of w, it is |lAfc(c— = |lA^c(c— < 
£ and therefore ||Mc(w - v)\\ < \\Mc{w - c)|| + ||Mc(c - < 2e. Likewise, 
it is straightforward to show that |lMc(c — c')|j < |jMc(c — v)\\ < e imphes 
||Afc(?7T- — c')ll — ^- Therefore, it is 



^2 ll^' - c 



||c — c 

.-2 ll^'-cl 



'^{w- vfMl [M^\Qc' ~ Qc)M~'] M,{m - c')| 



< Is^'f 'j^\\Mciw~v)\\p{M-\Qc' -Qc)M-') ||Af,(m-c') 

1 1 c — r' 1 1 

11?) — r.W 

< 2, 



!L_|p(M-*(Q,,-Q,)Af-i) 

Letting A = Mc'M'^ - I, it is M~\Qc> - Qc)M~^ = A* A + A + A^, and so 
p{M-\Q^, - Qc)M-^) = p{AY + 2p{A). By Lemnia[2) it is 

p{A) - p {{M,, - Mc)M-i) < a||A/,(c - c')|| 

and, using || Afc(c — c')|| < ||Afc(c — w)|| < e, it is 

- p^['^'l|A^c(c-c')ir + 2a||A/,(c-c')||] 

< [a2||A.4(c - + 2ct||A4(c - 

< {ecjf + 2ea 



which in turn imphes a <2 (ecr)^ + 4e(7. 



□ 



Lemma 4. Given an asymmetric e-net V , if v^w € V are Voronoi-neighbors, 

DW 

and it is q ^ i?y and c as described above, then it is 

(i = |£-2(c - v)'Q,{w -v)\>\l (2(1 + eaf) 

Proof. We first show that (w — vYQdv + w — 2c) ~ 0. Given that Dq {v, c) = 
Dq (w, c) imphes {v — cYQdv — c) = {w — cYQdw — c), it is 

{w-vYQciv + w-2c) = [{w - c) + {c~ v)Y Qcliw ~ c) + {v - c)] 

= [{w - cYQc{w ~c)-{v- cYQc{v - C)] 

+ [{w - cYQdv - c) + (c - w)*gc(«; - c)] = 

By Lemma |8] it is 

= |2(w - v)Qdc -v) + {w- vYQdv + w - 2c)| 
= 2\{w^vYQdc-v)\ 
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and thus \e-^{c-vYQc{w -v)\ > 1/ (2(1 + ecr)^). □ 

The next theorem uses the bounds of Lemmas [3] and [4] to prove that, under 
certain circumstances, the difference (f){c) — 4>{c') cannot vanish and therefore 
the anisotropic Voronoi diagram is orphan-free. 



Theorem 1. Given a continuous metric Q, the Du/Wang diagram of an asym- 
metric e-net (with respect to Dq ) is orphan free if ea < 0.09868. 

Proof. Given the construction at the beginning of this section, it must be (j){c) — 
(j){c') = 0. However, if ea < 0.09868, it is 



10(c) - <?!.(c')| > /3 - a > 1/ (2(1 + eaf) - 2 {eaY + iea > 

in^ c'll 

reaching a contradiction. 

Since ah points in qv must be closer to v than to any w ^ v, q cannot be 

DW 

in an orphan region of v. Additionahy, since for every point q S , the 
segment qv is also in , every Voronoi region is star-shaped with respect to 
its generating site. □ 

Finally, note that connectedness of Voronoi regions is shown by proving 
the stronger condition of star-shapedeness, which suggests that the condition 
ea < 0.09868 may be conservative in some cases. 

4.2 Orphan-free Labelle/Shewchuk diagrams 

In [S] , a condition is described for a two-dimensional LS diagram to be orphan- 
free. Although this condition is somewhat technical, an accompanying iterative- 
insertion algorithm is provided that, given enough time, will output an orphan- 
free LS diagram. In this section we describe conditions for an LS diagram to 
be orphan-free in any number of dimensions. The conditions are very similar to 



those of Section 4.1 namely, that the set of generating sites form an asymmetric 
e-net with respect to Dq , with sufficiently small e. The net requirement is 
somewhat natural in the sense that it implies that the sites are "uniformly 
distributed" with respect to Dq . 

Similarly as in Section |4.1[ we consider a point q that is in an orphan region 
of some V V. Since q is in an orphan region of f , the segment Tjv cannot be 
contained in , and thus there must be p' € that belongs to the Voronoi 
region of some w ^ v. This in turn implies that there are two distinct points 
c G qp' , c' G p'v that are equidistant from v, w. If we define the function 



LS 



mi)) = [Dq iw,p{l)) - Dq {v,p{l)) 
with p{l) = v{l — I) + cl, then it must be (f){c) = 4'{c') = 
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We now prove that q cannot be in an orphan region by showing that 4>{c) — 
0(c') cannot vanish, reaching a contradiction. 
It is 

= ^i^|0(c)-0(c')|=£-'^^|(c-cTQ«,(c'-^«)-(c-c')*Q.(c' 
1 1 c — c' 1 1 1 1 c — c' I 

= |£-2(c _ vfQ^{v ~w)+ £-2(c - vY{Q.u, - Q,){c' - v)\ 

which, letting a — \e^^ic — vY{Qu, — Qv){c' — v)\, and /3(Z) — \e~'^{c—vYQ.^{v — 
'w)\, can be rewritten, similarly as in Section 



4.1 



^^|<^(c)-</.(c')|>/3-a 
\\c-rJ\\ 

We now prove upper, and lower bounds for a and /3, respectively. 



Lemma 5. Given an asymmetric e-net V, ifv,w(£V are Voronoi-neighbors, 
and q G , and p{l) is as described above, then it is 

a ^ - vYiQ^ - Q,){c' - v)| < 7' + 27 

where 7 = £cr(l + k), and k = {1 + ecr) /(I — ea). 

it is 



Proof. Since ||Af„(c' — f)|| < ||Af^(c — < e, by Lemma 11 



a = \s'^{c - vYmIM-\Q^ - Q.,)M~HUc' - v)\ 

< e-^M,{c - 1^)11 ||M,(c' - v)\\p{M~\Q^ - Q,)M-^) 

< e-^M,{c-v)\\^p{M-\Q^~Q,)M-^) 

< (ecr(l + /c)) V 2ecr(l + fc) 7^ + 27 



□ 



Lemma 6. Given an asymmetric e-net V, ifv,w€V are V oronoi-neiqhbors, 

,1-3 , , , , , , „ „ I , J 

then it IS 



q £ Ry , and p{l) is as described at the beginning of Sec. 4-^ 

P = |e-2(c - vYQUy -w)\>{k^-j^- 27)/2 
where 7 £(7(1 + k), and k = {1 + £<t)/(1 — £cr) 
Proof. Because {v,c) = {w,c), by Lemma |llj it is 

\{w-vYQw{v + w-2c)\ = \{w - cYQw{w - c) - {v - cYQviv - c) 

- {v - cYiQw ~ Qv){v - C)\ 

= \{v~cY{Qn.-Q.){v-c)\<J^+2j 
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By Lemma [T0| it is 

e^/k^ < \{v-wYQ^{v-w)\ = \2{w~vfQ^{c-v) + {w-vYQ^{v + w-2c)\ 

< 2\{w- vYQ^,{c -v)\ + \{w- vYQ^{v + w-2c)\ 

< 2|(u;-w)*Q„(c-w)|+ 72 + 27 

and therefore /3 ~ \£^^ic — vYQwiv — w)\ > (fc^ — 7^ — 2^)/2, as elaimed. □ 
These bounds on a and /3 imply the following: 

Theorem 2. Given a continuous metric Q, the Labelle/Shewchuk diagram of 
an asymmetric e-net (with respect to Dq ) is orphan free if sa < 0.0584. 

Proof. Given the construction at the beginning of this section, it must be 4>{c) — 
(j){c') = 0. However, if ecr < 0.0584, letting 7 = ea{l + k), and k = {l + ea)/{l- 
£(t), it is 

fc^|0(c) ~ cf,{c')\ >l3~a>{k^-j^- 27)/2 - 7' - 27 > 

reaching a contradiction. 

Since all points in qv must be closer to v than to any w ^ v, q cannot 
be in an orphan region of v. Additionally, since for every point q € Ry , the 
segment qv is also in i?„ , every Voronoi region is star-shaped with respect to 
its generating site. □ 



5 Conclusion 

This paper presents a simple and natural condition for the two definitions of 
anisotropic diagrams of [S] and [S] to be composed of connected (and in par- 
ticular star-shaped) regions. The condition is simply that the generating sites 
be roughly "uniformly" distributed (forming an asymmetric £-net), with re- 
spect to the underlying metric. Apart from being natural, this condition is 
also commonly employed for certain practical problems. In particular, for opti- 
mal quantization, where we are interested in the primal Voronoi diagram, the 
optimal quantization sets form an e-net [7J [3j . For C°° PL approximation of 
functions [3], where we are interested in the dual simplicial complex, existing 
asymptotically-optimal constructions use vertex sets that form an e-net . 

Note that, although the definition of e-net must be slightly modified for 
the problem that concerns us here, this modification is not of great practical 
importance since existing algorithms for computing e-nets [B] are easily adapted 
to produce the desired (asymmetric) e-net. 

As mentioned in Sec. [sj computing an e-nct using the algorithm of [6 in- 
volves, at each iteration, computing the farthest point from the current set of 
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sites: a task equivalent in cost to computing the Voronoi diagram of each in- 
termediate set of sites. This makes exphcitly computing a unnecessary since 
computing an asymmetric £-net of sufficiently small e has a similar cost to 
iteratively running the algorithm of [5] while checking whether intermediate di- 
agrams are orphan-free, and then stopping when the first orphan-free diagram is 
produced. As mentioned earlier, the proofs in this paper guarantee that such an 
iterative algorithm will stop, while the specific bounds may give some indication 
of when this happens. 

In the eventuality that ways to compute (asymmetric) e-nets arise that are 
more efficient, it may be that computing a becomes useful since, along with 
Theorems [T] and [2] it would provide a simple lower bound on the largest e for 
which an asymmetric e-net is guaranteed to result in an orphan-free diagram. 
In this case it would become useful to know of more efficient ways to compute a. 
In particular, in [1] [this reference is to a supplementary document included in 
the submission, which will be cited as a Technical Report in the final version], 
we show that if the metric Q is continuously differentiable, more efficient ways 
to bound <T exist, which involve looking at every point in the domain only once, 
as opposed to comparing all pairs of points (as in Def. [2]). In particular, if the 
metric is specified as a PL function over a simplicial complex, then bounding a 
only requires computing a single number at each element (in constant time), and 
then taking the maximum over all elements, and has therefore linear complexity 
in the number of elements in the complex. 

Finally, note that, aside from the well-behaved-ness implied by the lack of 
orphan regions in the Voronoi diagrams, it is possible that orphan-freeness may 
be useful in guaranteeing properties of their dual abstract simplicial complexes 
such as being absent of inverted elements. In future work, we would like to 
explore whether orphan-freeness in a Voronoi diagram can be used, possibly 
along with additional conditions, to guarantee that it's dual is an embedded 
simplicial complex. 
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Appendix 

Assume given e > and a metric Q G C'^ defined over domain 0, and let 
k = {1 + ea)/{l — ea) > 1. The following lemmas are used in the proofs of 
Section m 

Lemma 7. Given a non-singular matrix A E M"^", it is p{A^^) — pm{A)^^ . 
Proof, li li, i — 1, . . . , n are the eigenvalues of A then 

p{A-') = max|/ri| = max|/,ri = {min\k\)-^ = p„M)~' 

i i i 

□ 

DW 

Lemma 8. Let V be an asymmetric e-net w.r.t. Dq , and v,w he Voronoi 
neighbors of the resulting DW-diagram. If c E is in the Voronoi regions of 
V, w then 

\\M,{v~w)\\ >£/(1 + £(t) 

Proof. Since V is an asymmetric e-net, it is either Dq {v, w) > e oi Dq {w, v) > 
e. Assume w.l.o.g. that Dq^ (w, w) > s and therefore |jM^,(i; — w)|| > s. Since 
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c IS m the Voronoi regions or v,w, it must be Uq (v,c) = Uq [w,c) < e. 
Therefore, by Lemmas [2] and [7j it is 

\\M,{v ~w)\\ = \\M,M-'M,iv -w)\\> p^{M,M~')\\M,{v - w)\\ 
= p{M^M-^y^\\M^{v - w)\\ >e/{l + e<j) 

□ 



Lemma 9. Let V be an asymmetric e-net w.r.t. Dq (resp. Dq ), and v,w G 
V be Voronoi neighbors of the resulting DW diagram (resp. LS diagram). Then 

l/k < p„(M^Af-i) < p{M^M-^) < k 

Proof. Since v,w are Voronoi neighbors, there is a point c G that befongs to 
the Voronoi regions of both v and of w. In an DW diagram, it is |j Mc(c — v)\\ = 
\\Mc{c — w)\\ < e, and therefore, by Lemmas [2] and [Tj it is 

p(M^M-i) < p{M^M-^)p{M,M-^) 

= p{M^M-^)p^(M^M-^)-^ < = k 

1 — ecr 

and 

PmiM^M-^) > Pra{M^M-^)pra{M,M-^) 

= p„,(Af^Af-i)p(M„M-i)-i > = l/k 

1 + £(T 

In an LS diagram, it is ||Mt,(c — v)\\ = \\Miu{c — w)\\ < e, and therefore, by 
Lemmas [2] and [T) it is 

p(M„M-i) < p(M^M-^)p{M,M-^) 

= p„(MeM-i)-V(McM-i) < = k 

and 

p^{M^M-^) > p,n{M^M~^)p^{M,M-^) 

= p{M,M-')-'p„,{M,M-') > 1^ = l/k 

1 + ea 

□ 



Lemma 10. LetV be an asymmetric e-net w.r.t. Dq , andv^w G V be Voronoi 
neighbors of the resulting LS diagram. Then 

e/k < \\My{v - w)\\ < e(l + fc) 
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Proof. Since v, w are Voronoi neighbors, there is c S that is in the LS Voronoi 
regions of v,w, and therefore satisfies ||-Mi,(w — c)|| = ||M^(w — c)|| < e. Thus, 
by Lemmas [2] and |9) it is 

\\M,{v-w)\\ < \\AMv-c)\\ + \\M,{c-w)\\ = UMv-c)\\ + \\M,M;^^ M^{c ~ w)\\ 

< WM^v - c)\\ + p{M,M-^)\\M^{w - c)\\ < e {1 + p{M,M-^)) < e{l + k) 

and 

UUv - w)\\ > \\M,M-'M^iv ~ w)\\ > PmiM,M-')s > e/k 

□ 



Lemma 11. LetV be an asymmetric e -net w.r.t. Dq , andv,w€ V be Voronoi 
neighbors of the resulting LS diagram. Then 

p[M-'Q^M-^ -I)< (£(7(1 + k)f + 2e<7{l + k) 

Proof. Let A = M^M'^ - / and B = M-*Q^M-^ - I, where it is B = 



it is 



A* A + A + A*. By the definition of a and Lemma 10 

piA) = p{M^M-^ -I)< n\\M,{v - w)\\ < ea(l + k) 
and therefore 

p{M-'Q^M-^-I) = p{A*A+A+A') < p{Af+2p{A) < {ea{l + k)f+2e(j{l+k) 

□ 
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